XVI International Symposium on Very High Energy Cosmic Ray Interactions 
ISVHECRI 2010, Batavia, IL, USA (28 June 2 July 2010) 



1 



Physics of high energy atmospheric muons 
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In the first part of the talk the interesting new results of MINOS, OPERA and CMS collaborations (connected 
with the observational evidence of the rise in the muon charge ratio at muon energies around 1 TeV) are 
briefly discussed from theoretical point of view. A short review of charge asymmetric effects in muon energy 
losses is given. In the second part of the talk, the modern theoretical approaches to the problem of heavy quark 
production in high energy nucleon- nucleus interactions are briefly considered (color dipole formalism, saturation 
models). The recent new theoretical developments in the ancient problem of intrinsic charm are also discussed. 
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1. THE CHARGE RATIO OF 
ATMOSPHERIC MUONS 

1.1. Pika model 

A simple formula for the muon atmospheric charge 
ratio is obtained from the approximate expression for 
the differential atmospheric muon spectrum [l|, : 



dE„ 



OiZNi 



dE^ ^ l + b,E^cose*/ei 



(1) 



where d4)N/dE^ = d(j)N/dE (at E ~ £'^t) is the pri- 
mary spectrum of nucleons (E^ is the muon energy in 
atmosphere), index i numerates muon parents (tt, K 
etc), the constants a^, hi are purely kinematical (they 
depend on the characteristics of parent's decay into 
muons), ei are the critical energies, which are defined 
as parent's energies in atmosphere, at which their in- 
teraction lengths are equal to decay lengths (for a ver- 
tical propagation), 9* is the zenith angle of the muon 
trajectory (for simplicity, we assume that atmosphere 
is flat). At last, Zij is the spectrum weighted moment 
defined by the expression 



Zij — 



1 daij 



(2) 



where CTy is the inclusive cross-section for the produc- 
tion of a particle j from the collision of a particle i 
with a nucleus in the atmosphere and xiab — Ej/Ei is 
the energy fraction carried by the secondary particle, 
7 is a spectral index of the primary nucleon spectrum, 
7 = 1.7. 

If we are interested in the interval of muon energies 
(1 — 10"*) GeV, one can take into account only pions 
and kaons as muon's parents. In this case, the muon 
spectrum can be expressed as 



RKiTaKfK± 



dEf, l + b^E^cos0*/e^ I + bxE^cose* / ex 



(3) 



where Rkt, = Znk/Zn-,t and = 1 - = 

Zntt+IZntt, fK+ = 1 - Jk- = Zpfii+fZNK are frac- 
tions of all pion and kaon decays that yield posi- 
tive muons. This simple formula used for an anal- 
ysis of atmospheric muon charge ratio data is a ba- 
sic of the so-called "pika model" More exactly, 
this model uses the following main assumptions: 1) 
particle ratios, 7r+/7r~ = and / K~ = fK+, 
are energy independent; K/n ratio is given by the 
Rktt ■ aK/a^ = 0.149 • 0.246/0.674 = 0.054 ^, 2) all 
the spectrum weighted moments, Zij, are energy inde- 
pendent, 3) the muon charge ratio does not depend on 
the Ef^ and the cos 9* separately, but on the product 
Ef^ cos 9*, where i?^ is the energy of vertical muons 
on the Earth surface, 4) contributions from charmed 
mesons are ignored. 



1.2. Particle ratios 



It is well known from studies of iiT-production in 
NN and A^yl-coUisions, in the energy region near 
the threshold, that K^/K^ ratio is rather large. 
The explanation is simple: K'^-meson contains s~ 
quark whereas, for a final baryon, s-quark is needed. 
The transfer of s-quark from X^-mcson to a final 
baryon is realized most effectively through the reac- 
tion K~N — )■ ttA. So, as a result, the escape of K~- 
mesons is hampered. Measurements of particle ratios 
at high energies, far from the threshold of K pro- 
duction, shows (see, e.g., the data of the BRAHMS 
Collaboration [J| for pp-coUisions at ^/s = 200 GeV) 
that, at large rapidities (close to the fragmentation 
region), the K ~^ / K~-t aiios are systematically larger 
than 7r"'"/7r~-ratios. In general, such asymmetries in 
the production of particles and antiparticles are indi- 
cations of the leading particle effects in hadronic col- 
lisions and are rather intensively studied now experi- 
mentally as well as theoretically. In particular, large 
values of /K^ -ratios measured by [3] are well ex- 
plained by the DPM JET model [1]. 
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1.3. Muon energy losses: 2;-odd effects 

1 .3.1 . Ionization energy losses 

As is known, the standard expression for energy loss 
contains, except of the main term, proportional to 
z^, the tiny term which is proportional to (a 
is the fine structure constant) [gl. This term is due to 
the next-to-leading order correction to the relativistic 
Rutherford formula. The z'^-correction to dE/dx is 
taken into account in an analysis of data in the MI- 
NOS experiment |3| as well as in the CMS , in spite 
of its smallness. For the energy loss in matter is 
about ^ 0.15% higher than for for CMS experi- 
ment, and it affects the muon charge ratio by less than 
~ 0.3% over the entire energy range Q- 

1 .3.2. Radiative energy losses 

The muon brcmsstrahlung cross section is given by 
the formula 



da 



da 
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Here, uj is the photon energy, the first term in r.h.s. 
is the cross section calculated in the Born approxi- 
mation (including, in this approximation, the nuclear 
formfactor correction), Ac is the "non-Coulomb" cor- 
rection Q (it is important only for muons and is the 
non-Born part of the correction which arises due to the 
non-Coulombic character of the electromagnetic field 
of the nucleus having the extended size), Adbm is 
the well-known Coulomb correction of Davies, Bethe 
and Maximon jflj. At last, A.nonieading takes into ac- 
count the main correction to ultrarelativistic approx- 
imation which had been used before for a calculation 
of all components of the bremsstrahlung cross section. 
It appears [lOj that this correction is proportional to 
— -pT^Za, the minus sign corresponds to the u~ , i.e., 

the /i+ has again a slightly larger energy loss. So, the 
difference between the fi'^ and fi~ energy losses due to 
bremsstrahlung is inversely proportional to the muon 
7-factor, 
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and is very small for the case of the rock {Z — 11). 
This difference is negligible Q in the region around 
E^ ^ (0.1 — 1) TeV, where radiative energy losses 
become comparable with ionization energy losses. 

One should note that, of course, there are other 
sources of the charge asymmetry in muon energy 
losses. For example, in the process of e~''e~-pair pro- 
duction by muon on atomic electrons such an asym- 
metry arises due to interferences of different diagrams 
contributing to the process [lH. One should men- 
tion, in this connection, also the work fl^ in which 
diffractive corrections to muon bremsstrahlung had 
been considered. 
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Figure 1: R^, values measured by OPERA, MINOS and 
other experiments (taken from jl^]). 



1.4. Conclusions 

1. The recent experimental results of three groups, 
MINOS [|, OPERA [11] and CMSQ (together with 
the old result of Utah detector [1J|) show that the 
atmospheric muon charge ratio at energies around 1 
TeV and higher is sizeably larger (i?^ ^ 1.36 — 1.4) 
than at energies around ~ 10^ GeV (see Fig. [Ij. 

2. The data fittings give the similar results for all 
three groups: r^ = - U+) = 1.241 ± 0.035, 

= fK+/{l-fK+) = 2.26±0.29 (MINOS N-hFjl); 
= 1.229 ± 0.001, = 2.12 ± 0.03 (OPERA^|); 
U+ ^ 0.553 ± 0.005, fK+ = 0.66 ± 0.06 (CMS 



The increase of the ratio i?,, 



is due to the 
fact that i) according to the pika model, at ener- 
gies ^ 1 TeV, the contribution of kaons in i?^ be- 
comes dominant, and ii) for kaons, the ratio K'^ / K~ 
is larger than the corresponding ratio for pions. 

4. The fits of the pika formula to the old results of 
calculations [ll| and simulations [H, lead to the 
following numbers 0: = 2.87, r„ = 1.25 (for the 
calculation by Lipari [lal); fK = 2.39, = 1.26 (for 
the simulations with CORT code [la]); r^ = 1.63, 
r^r = 1.28 (for the simulations by Honda [17|). 

5. It would be quite interesting to check whether 
the recent experimental results [3| and theoretical pre- 
dictions [^ are consistent with the new i?^-data. 



2. CHARM PRODUCTION IN 
pA-COLLISIONS AND PROMPT MUONS 

2.1. Introduction 

One can see from Fig. [T] that the value of Rf^ at 
Efi ^ 10 TeV strongly depends on the charm pro- 
duction model. It is well known that it is not so 
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a) 



b) 




c) 




Figure 2: The schematic diagrams of the pyl-collision. 



difficult to calculate muon atmospheric fluxes, but 
only if you know well the fluxes of niuon parents 
(charmed mesons, in a given case). So, the physics 
of hadronic charm production is intimately connected 
with physics of atmospheric high energy muons. 

The schematic diagram of the process is shown in 
Fig [2^. We consider in this talk only two examples of 
such calculations. 

1. If one uses perturbative QCD approach, the typ- 
ical diagram for the lower blob of diagram a) is shown 
in Fig. [2b (gluon fusion). 

2. If it is assumed that there are non-perturbative 
"intrinsic" heavy quark components in the proton 
wave function 18], the diagram in the lower blob of a) 
describes the inclusive scattering of the heavy quark 
of the proton projectile on the target (Fig. [5J;). 

The common points in these cases are following: 
the fragmentation region of the proton projectile is of 
interest; the target is a nucleus; the proton energy is 
very high ( > lOOTe^). We will show in this talk that 
both cases can be studied using the modern approach 
to high energy scattering in QCD based on the idea 
of the Color Glass Condensate (CGC) (for a review, 
see [H). 

2.2. CGC framework 

2.2.1 . Dipole cross section 

In the CGC approach, the proton and the target nu- 
cleus are effectively described as static random sources 
of color charge on the light cone. The leading contri- 
bution to particle production in pA-collisions is ob- 
tained by calculating the classical color field created 



by these sources and then averaging over a random 
distribution [20*1. 

The field is a solution of the classical Yang-Mills 
equation: 



(6) 



where J'^ is the classical current which is determined 
by the sources Pp and p^: 

j: = g5''+5{x-)pp,a(:xi_) + g5''-5{x+)pAA^^)- (7) 

Here, — ^ ^ , x± is the transverse component of 

X, x^ is the "LC time". Pp,pA are static, i.e., inde- 
pendent of the LC time. S{x~) function assumes an 
infinitely thin sheet of color charge. 

It is shown in many works that in CGC approach 
(see, e.g., [2l|) the amplitude for a scattering of a 
quark from the target, qA qX, is expressed through 
the "gauge rotation" matrix U{x±): 



MqA^qX = u{q) X 

This matrix is given by the expression: 



u{p). (8) 



U{x±) — Texp 



dx' 



-^PA,a{x^,X±)t°' 



(9) 

is a color matrix, generator of the fundamental rep- 
resentation of SU{Nc). Here, pA,a{x) is the density of 
color sources in the nucleus in the covariant gauge 
dfj.A^ = 0. For concrete calculations the LC-gauge, 
in which ^+ = 0, is used, because in this gauge the 
color fields depend only on the matrix U (i.e., they are 
"pure gauge" fields). But U depends on the pA,a{x), 
the density of sources in the covariant gauge. 

It is very important that the same matrix U{x\^ 
enters the expression for the dipole-hadron cross sec- 
tion, crdip(^_L)7 which is the total cross section of in- 
teraction with nucleus, for a dipole of transverse size 
r*! = x^_ — y± (with the quark located at xj_ and the 
antiquark at y± ) and is written as an integral over 
all the impact parameters X± = {x± + y±)/2, 



o-diplr^ 



-u(x^_ 



(fXiTr^{l 

U+IX, 



(10) 



We will see, in the next Section, that charm pro- 
duction cross sections contain just the crdip(^_L)i 
which is expressed, as is seen from Eq. pUj) . 
through the 2-point function, or "2-point correlator" , 
{U{x±_)U+{yj_)). In the Eq. Tr^ is the color 

trace {U is the matrix in the color space), and en- 
semble average is performed over the sources. In the 
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Figure 3: The forward dipole amplitude, (Tdip/27ri?^, as a 
function of the transverse separation between g and g 
(taken from :22i]). 



McLerran-Venugopalan (M-V) [201 model a Gaussian 
distribution of color sources is used, 



W\p\ — exp 



(O) 



2p\{x+) 
[dp]W[p]0[p]. 



(11) 
(12) 



The color source is "frozen" through the collision due 
to time dilation but changes from the collision to the 
subsequent one. Evidently, p\(x'^) is a density of 
color sources per unit volume and 



(13) 



y-A = / dx-^pj^ix-^) 



is the average density of the sources per unit trans- 
verse area, 



{Paix±)pb{y±)) = Sab5{x± - y±)p\, 



2 

l^A 



. (14) 

The straightforward calculation gives for adip{r±) the 
result (if the target is homogeneous in the transverse 
plane) : 

adip(rl) - 2^i?2 _ (_g2/(^^))] , (15) 
where f{r±) is the known function and is given by 

Q^ = ^(tX)Mi = ^4^M^-a?^^/^ (16) 

By definition, is the saturation scale (in the M-V 
model). The schematic form of crdip(^_L) is shown in 

Fig. m 

Taking into account the quantum corrections |19j . 
one can show that the dipole cross section can be pa- 
rameterized in the following form [23j : 

adip(r^,:r) « 2nR\ (^1 - eM-^-^^^^^)) , 

(17) 



1 GeV^ 



2N, 



xo « 3 X 10"*, A « 0.3, 7 « 0.628. 



(18) 
(19) 



It follows from this equation that the saturation scale 
Ql becomes dependent on energy, because the dipole 
cross section depend now on Bjorken x. In DIS, the 
value of X is connected with the photon virtuality, 



2Mqo 



(20) 



In a case of the inclusive hadronic scattering and, in 
particular, in the case of heavy quark production in 
hadronic collisions, x is the momentum fraction of the 
target parton participating in the reaction (x — xa)- 

2.2.2. Kinematics of reactions 

1. In the case of intrinsic charm production (Fig. 
[5J:), there is the simple connection between xa and 
the momentum fraction of the projectile (xp): 



XA 



Xa 



(21) 



where aij? is the Feynman x of the produced particle; 
yh is its rapidity, qt is its transverse momentum. In 
the forward region one has 



a^i?— J-l; i < Xp < xf, and 



{xa) -^^^0- 

(22) 

2. In the case of extrinsic charm production (Fig. 
[21d) we have two colliding gluons, with momentum 
fractions Xg and xa- Assume, for simplicity, that the 
heavy quark pair is produced with Feynman xp cor- 
responding to a fixed Mqq. One has 



XA 




XF + XA, 



(23) 

and, again, in the forward region, xa ^ ^'^qqI ^ ^ ^■ 

2.2.3. Target description in a small x region 

1. Coherent lengths of small x partons in the tar- 
get are very large, i ~ 1/{Mmx), much larger than 
the nuclear size. Therefore, at very small Bjorken x 
inclusive particles are produced in hadron-nucleus col- 
lisions coherently by the nuclear color field. 

2. The gluon distribution which is the number of 
gluons per unit rapidity in the hadron wave function. 



xG{x,Q'^) = dNglnons/dy, 



(24) 



grows with the decrease of result, the satura- 

tion scale Ql also grows: 



(25) 
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Physically, is the average color charge squared of 
the gluons per unit transverse area per unit rapidity. 

3. If the density of color charge grows, the classical 
description of the color field becomes possible as a 
first approximation [20l |. The nonlinear effects in the 
classical Yang-Mills equations and quantum evolution 
corrections lead to a saturation of gluon distributions. 

One should note that the experimental as well as 
theoretical study of heavy quark production at for- 
ward rapidities is very important just for the "satura- 
tion physics" because the saturation scale Qs increases 
with energy, Qs ^ A^/^s'^/'*, and, at large A and s, 
can be larger than Mq. The typical transverse mo- 
menta of partons in the hadron are of the same order 
as Qs , so, it is natural to expect that in such cases 
the probability of QQ-pair production will be greatly 
enhanced. 

2.3. Extrinsic and intrinsic charm 
production 

2.3.1 . Extrinsic charm in tlie color dipole approach 

In the model of [l^, HH] the differential cross section 

is 

^ (26) 

(27) 

Here, |^'g(-2^i is the probability of finding QQ- 

pair with a separation rj_ and a fractional momentum 
z in the gluon, /i ~ mq is the factorization scale. The 
cross section is written in a factorized form and con- 
tains (JdipixAiri^ which depends only on two-point 



correlator {U{x±^U^{y±y). It is a consequence of the 
special choice of perturbative diagrams. In general, if 
the nuclear color field is strong, one must take into 
account all orders in pA (if g^PA ^ 1) and, in par- 
ticular, include diagrams with exchange of two gluons 
[26J. After this there will be no simple factorization: 
the cross section will depend on 2-point, 3-point and 
4-point correlators. 

On Fig. m the predictions of this model (131 are 
shown, for several incident proton energies, with a 
comparison with the corresponding results of pQCD 
model [i^. One can clearly see the suppression of 
charm production at very high energies due to gluon 
saturation effects. On the Fig. [SJ the prediction of 
the same model [13] for the prompt lepton fluxes are 
shown (in comparison with the conventional fluxes 

Hi). 

2.3.2. Intrinsic charm production. 

In the CGC approach, the diagram of the inclusive 
heavy quark scattering (Fig. [5J:) is the basic one. The 
amplitude of this scattering is given by Eq. 

After squaring this amplitude, one obtains, again, 
the formula containing 2-point correlator, which can 
be expressed through the dipole-hadron cross section 
(Tdip- The formula for the inclusive production of the 
charmed meson is [13, HH : 

7 n ~ 77, Y? ( dXp fc/pi^piQ ) X 

dyd^p±_ (27r)^ J xp ' 

X / d\^e^^-^-- "''^ 7 ' ""^^ DnU^, Q') . (28) 

The existence of an "intrinsic charm" (IC) compo- 
nent of the nucleon is not forbidden. Many nonpertur- 
bative models predict such a component at an energy 
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Figure 4: The comparison of the predictions of the dipole 
model [531 with the corresponding predictions of the 
pQCD model [2^, for different proton energies (taken 
from [131). 



10" 

7 10" 

'e 



10" 



*io- 



10" 









^■•^ Prnmpt 











4 5 6 7 

log,oe [GeV] 

Figure 5: The prompt muon flux predicted by dipole 
model [231 {Vfj. + i>i_i, Ve + Ve and fi^ -\- p~ - fluxes are 
approximately equal). Conventional fluxes are 
represented by dashed lines [2§]. The figure is taken 
from [131. 
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Figure 6: The comparison of a;F-spectra for BHPS (solid) 
and sea-like (dashed) models, for two energies, -y/s — 10^ 
and 10'' GeV {A = 14, from g^). 




Figure 7: The comparison of the results of [34I (thick) 
with the results of [3J| (thin), for Ep = 10® GeV (uppe 
curves) and 3.4 x 10* GeV (lower curves); A = 14. 



scale comparable to mass of charm quark m^. In |32| . 
the predictions of the specific light-cone model devel- 
oped in [1^ ("BHPS model") are used for the calcu- 
lation of charm production (authors of [l^ calculated 
the PDF of charm quarks from the 5-quark component 
uudcc in the proton). As a second example, authors 
of [3?! considered the phenomenological model [ssj . 
in which the shape of the charm PDF is sea-like, i.e., 
similar to that of the light flavor sea quarks (except 
for normalization). In the concrete calculation, the 
scale dependence of PDFs must be taken into account. 
For this aim, authors of 32] used the PDFs from the 
work [33], based on the calculations of CTEQ group. 
For BHPS model, it was assumed, in js^, that the 
charm content of the proton is on the maximal level, 
{x)c+c = 0.02. The same is for the sea-like IC model: 
(x)e+e = 0.024. 

On the following two figures some results of the re- 
cent work {3^ are presented. Fig. [5] gives the com- 
parison of X i?-spectra for two different models of the 
intrinsic charm. On Fig. [7] the comparison of the pre- 
dictions for inclusive cross sections for intrinsic [3^ 
(the total contribution of all mesons) and extrinsic 
charm 34 1 is shown. Authors of [s^] used the dipole 
model approach and their results are similar with the 
resuhs of pTj . 

It follows from figure [7] that in the tail region, xp > 
0.7, the curves of [S^l are even steeper than those of 
[34| , in spite of their "intrinsic" origin. One must bear 
in mind, however, that in [s^l the inclusive spectra of 
charmed hadrons rather than cc pairs are calculated. 
Naturally, a taking into account of the fragmentation 
of c quarks into hadrons leads to additional steepening 
of XF-spectra. 



2.3.3. Experimental data of RHIC 

One can see from the Fig. |8] that the results of 
different pQCD approaches (presented on the Figure 
by curves) underestimate, in general, the cross section 
of charm production at high energies. Authors of [ssl ] 
argue, that it may be connected, in particular, with 
large non-perturbative contributions in high density 



- P(p) + P cc 












' 1 1 II 




0' 
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Figure 8: The total cross section of charm production in 
pp collisions at p > 400GeV/c [s^. At high energies, 
squares correspond to PHENIX data and triangles to 
STAR data. The data of cosmic rays are shown by 
diamonds. The figure is taken from [35i] . 
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states (e.g., with CGC effects). 

2.4. Conclusions 

We tried to show in this review that the raodern 
approach to studies of proton-nucleus nucleus inter- 
actions based on the CGC theory can be effectively 
used in calculations of heavy quark production in pA 
collisions. The most important application of this ap- 
proach is an use of it in calculations of open charm 
production in collisions of cosmic rays in atmosphere 
(needed for subsequent calculations of prompt muon 
and neutrino fluxes at energies around 100 TeV and 
higher). The reason is simple: in collisions of cosmic 
rays in atmosphere, the fragmentation region of the 
projectile (or forward rapidity region) is most impor- 
tant for calculations of all spectra of secondary parti- 
cles, and just in this region the effects of gluon satu- 
ration studied in CGC theory are most efficient. 
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